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16.  Abstract 

— ^To  obtain  the  aircraft  ground  speed  during  Instrument  Landing  System  approaches  an 
algorithm  is  developed  which  utilizes  the  glide  slope  deviation  signal  and  the  rate  of 
descent  as  determined  from  the  radar  or  barometric  altimeter.  The  accuracy  of  these 
inputs  is  determined  from  a study  of  all  appropriate  Federal  Aviation  Administration 
Flight  Inspection  Records  and  from  examination  of  topographic  maps  of  the  areas 
beneath  the  flights  for  which  records  were  available.  It  was  found  that  irregularities 
inlhe  glide  slope,  variations  in  the  slope  of  the  terrain,  and  rapid  fluctuations  in  ihe 
barometric  pressure  associated  with  wind  shear  will  all  induce  excessive  errors  in  the 
value  obtained  for  the  ground  speed. 
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INTRODUCTION 

An  aircraft  flying  at  low  altitude  and  airspeed  during  a landing  approach 
is  vulnerable  to  any  sudden  change#  in  wind  velocity,  which  at  best  makes  main- 

K 

taining  a controlled  descent  difficult,  and  at  worst  may  result  in  a loss  of  lift 
sufficient  to  cause  a crash.  It  has  been  determined  that  the  availability  of  airspeed 
and  ground  speed  to  the  pilot  is  a substantial  aid  in  managing  the  craft  under  these 
conditions  of  wind  shear.  ^ 

for  economic  reasons,  a method  of  ground  speed  determination  involving 
equipment  already  present  on  commercial  aircraft  was  found  to  be  a desirable 
objective.  To  this  end  an  investigation  was  made  into  the  possibility  of  using  the 
radar  or  barometric  altimeter  derived  descent  rate  plus  the  glide  slope  deviation 
signal  to  calculate  ground  speed  during  an  Instrument  Landing  System  (ILS)  approach. 
All  appropriate  FAA  Flight  Inspection  Data,  consisting  primarily  of  chart  recordings 
produced  by  the  Arma/FAA  Automatic  Flight  Inspection  System  were  obtained  from 
the  FAA  offices  in  Oklahoma  City.  Data  involving  the  barometric  altimeter  was 
unavailable  for  study,  but  meterological  conditions  associated  with  wind  shear  were 
found  to  impose  serious  doubts  as  to  the  value  during  wind  shear  encounters  of  any 
barometric  altimeter  based  ground  speed.  Also,  United  States  Geological  Survey 
Topographic  Maps  were  studied  to  assess  the  effects  of  slope  of  terrain. 


To  summarize  the  results  of  the  investigation,  the  proposed  method  is 
inadequate.  The  effects  of  glide  slope  irregularities  and  slope  of  terrain  are 
so  adverse  that  further  consideration  is  unwarranted.  It  is  therefore  recommended 
that  other  techniques  of  ground  speed  determination  be  identified,  studied  and 


evaluated. 


CHAPTER  II 


F 
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GROUND  SPEED  ALGORITHM 


A runway  equipped  with  an  instrument  landing  system  gives  vertical 
guidance  along  the  proper  descent  angle  to  the  correct  runway  "touchdown" 
point  by  means  of  a radio  beam  directed  from  the  glide  slope  antenna  outward 
along  the  approach  course  at  about  3 degrees  above  the  horizontal . The  beam  is 
modulated  at  two  different  frequencies,  the  equal  signal  region  defining  the 
proper  descent  path,  and  the  deviation  from  equality  indicating  the  angular  dis- 
placement of  the  aircraft  above  or  below  the  glide  slope.  A pilot  making  an 
instrument  landing  approach  can  determine  by  consulting  his  glide  slope  indicator 
the  angle  between  the  horizontal  plane  and  his  craft  as  measured  at  the  antenna. 

If  the  average  inclination  of  the  center  of  that  particular  glide  slope  is  given  as 
6 and  if  the  glide  slope  deviation  signal  is  6 then  the  elevation  angle  will  ideally  be 
0 + 6 . The  antenna  is  typically  located  about  1000  feet  from  the  approach  end  of 
the  runway  and  offset  400  to  600  foet  (d„)  to  one  side  of  the  runway  centerline. 

Thus  a pilot  approaching  directly  above  the  imaginary  extended  centerline  and  main- 
taining a constant  (normally  zero)  glide  slope  deviation  will  be  guided  toward  a 
touchdown  on  the  runway  abeam  the  glide  slope  antenna  (point  O in  Figure  1). 

Upon  approaching  or  crossing  the  runway  threshold  the  geometrical  effect  of  the 
antenna  offset  becomes  evident  and  the  use  of  6 for  guidance  must  be  abandoned. 

It  has  been  proposed  that  this  knowledge  of  the  aircraft's  position  with  respect 
to  the  glide  slope  combined  with  the  rate  of  descent  as  determined  by  the  radar 
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(or  barometric)  altimeter  be  used  to  calculate  the  ground  speed.  The  following 
relations  exist  among  the  quantities  illustrated  in  Figure  1: 


tan  4>  = (1  + (d0 /x)2)1^2  tan  (0+6)  = y/x 


2 2 

(x  +dQ)  =y/tan(0  + 6) 


It  must  be  emphasized  that  y is  the  altitude  of  the  craft  above  a plane  tangent  to  the 
Earth's  surface  at  the  glide  slope  antenna  and  x is  the  distance  to  the  craft  from  the 

n 1 / r\ 

runway  centerline  abeam  the  glide  slope  antenna.  The  quantity  (1  + (d0/xf)  c 
in  equation  1 is  not  precisely  one  because  the  angle  (J>  is  referenced  to  the  poi^  O, 
rather  than  the  base  of  the  antenna.  Differentiating  equation  2 with  respect  to  time 
and  d. noting  J • 


v = (1  + (d0/x)  ) (v  /tan  ( 0 + 6 ) — y 6 /sin  2 ( 0 + 6 ) ) 
x y 


If  dc  = 400  feet  the  multiplicative  factor  has  the  value  1 .02  for  x = 2000  feet, 

1 .005  for  x — 4000  feet,  and  1 .002  for  x = 6000  feet.  For  a landing  approach 
speed  of  250  ft/sec,  neglecting  this  correction  would  result  in  ground  speed  error* 
considerably  greater  than  5 ft/sec  in  the  final  1000  feet  to  the  runway  threshold. 

For  the  usual  antenna  array  the  glide  slope  is  not  a vertically  oriented 
cone  with  a half-angle  of  90  -0  and  vertex  at  the  antenna  base,  but  rather  it  is 
a hyperboloid  of  revolution  having  this  cone  as  its  asymptote  and  typically  having 
a 1 .5  foot  closest  approach  to  the  antenna  base.  To  allow  for  this,  the 


square  root  in  equation  3 should  contain  the  additional  term: 


( — — - )«(  -^-)2 
v x tan  6 ' v x 


For  x = 1000,  that  is  for  the  aircraft  not  yet  to  threshold,  this  changes  the  square 
root  by  less  than  .0005  and  this  correction  may  be  safely  ignored. 

Since  v^  will  be  obtained  from  either  radar  or  barometric  altitude,  which 
are  referenced  to  the  Earth's  surface  and  to  mean  sea  level  respectively,  correction 
for  the  curvature  of  the  Earth  is  required.  For  the  distances  concerned  the  following 
cpproximation  is  sufficiently  accurate:  (See  Figure  2) 

(4)  C = .88  C2 

C is  the  correction  in  feet  to  be  subtracted  from  the  radar  or  barometric  altitude, 
and  D is  the  distance  in  nautical  miles  from  the  glide  slope  antenna  to  the  craft. 

The  error  induced  in  the  vertical  velocity  is  found  by  differentiating  the  above  with 
respect  to  time: 


3f  - L76  C JT  ='i-77XI0-8xvx 


The  derivative  is  the  ground  speed  in  nautical  miles  per  second  and  the  constant 
in  the  second  half  of  the  expression  arises  from  conversion  to  feet  and  feet  per  second. 
For  a ground  speed  of  250  ft/sec  and  C = 4 nautical  miles,  the  correction  is  .29  ft/sec 
corresponding  to  a ground  speed  error  of  about  5.5  ft/sec  . For  D = I nautical  mile  it 
is  .07  ft/sec  corresponding  to  1.4  ft/sec  in  v . Both  this  and  the  factor  involving 


XLW-  L . - , ...  . 


dG  in  equation  3 require  knowledge  of  x,  which  is  available  from  the  Distance 
Measuring  Equipment  (DME).  Even  though  the  DME  is  not  sufficiently  accurate  to 
derive  a useful  ground  speed,  using  it  for  position  in  these  relatively  small  corrections 
gives  adequate  results.  For  example,  at  x = 2000  feet  or  about  1000  feet  from 
threshold,  assuming  x is  correct  to  within  .02  nautical  miles,  the  error  in: 


(1  +(d0/x)V/2 


is  limited  to  .3%  or  less.  Thus,  several  complications  which  cause  relatively  minor 
systematic  errors  have  been  noted  and  corrected.  The  complete  equation  for  v^  is: 


i / l Z7T7  4.77  X 10 

vx=  1/V'  + +T5TTTT T. 


V 

( tan  1 0 + 


The  data  available  for  analysis  is  the  output  of  the  Arma  Flight  Inspection 

System  for  a series  of  ILS-3  runs  on  three  different  glide  slopes,  RGR  and  OKC  in 

3 

Oklahoma  City,  and  PVO  at  Atlantic  City  (See  Figure  3).  In  these  runs,  the 
test  aircraft  flew  down  the  glide  path  producing  time  correlated  trace  recordings 
of  the  glide  slope  deviation  signal,  the  localizer  signal  and  of  the  altitude  as 
determined  by  an  inertial  sensing  system  and  also  as  determined  by  trie  radar  altimeter. 
Further,  for  every  .1  nautical  mile  traveled,  as  determined  by  the  INS,  an  indicator 
mark  was  placed  on  the  trace  (Figure  4).  The  usefulness  of  these  marks  is  limited 
to  giving  the  approximate  distance  from  threshold  to  aircraft  and  a rough  idea  of  the 
ground  speed  for  several  reasons: 

1)  The  marks  are  produced  for  every  .1  nautical  mile  traveled  or  only 
about  every  2.5  seconds. 


I 


CHAPTER  III 


GLIDE  SLOPE  IRREGULARITIES 


Since  v^  is  derived  from  y and  6 is  derived  from  6/  if  i*  in  general  necessary 


when  calculating  the  probable  error  in  v to  consider  the  correlation  between  the 


error  in  y and  the  error  in  v and  the  same  for 6 and  6 . This  results  in  a complicated 

/ 


expression,  which  only  obscures  the  arguments  made.  Therefore,  when  altitude 
data  is  considered,  a perfect  descent  is  assumed,  that  is  6=  0,  and  the  second  term  in 
equation  3 is  neglected  yielding  the  following  expression  for  the  probable  error  in  v^: 


(6) 


PE  (v  ) = PE  (v  )/tan  ( 0 + 6 ) 
x y 


This  is  reasonable.  If  the  altitude  data  is  not  sufficiently  accurate  for  other  wise 
ideal  conditions  it  will  never  be  adequate.  Similarly,  when  considering  the 
error  in  6 and  6,  it  is  obvious  that  an  error  in  6 has  relatively  little  impact  on 


v^,  because  normal Iy6is  much  smaller  than  0 and  occurs  only  in  the  combination  0 + 6 . 


To  a good  approximation: 


(7) 


PE  (v  ) - PE(6)y/sin  (0+6) 


Thus  attention  is  directed  primarily  to  the  errors  in  6 and  v^,  and  it  turns  out  that 


the  most  essential  problem  is  not  the  accuracy  of  the  altimeter  or  glide  slope 
receiver,  but  the  quality  of  the  glide  slope  and  the  nature  of  the  terrain  flown 


over. 
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According  to  the  U.S.  Standard  Flight  Inspection  Manual  the  glide  slope 
elevation  or  alignment  ( 0 ) may  be  permitted  to  vary  as  much  as  + .22°  or  + 7.5% 
from  its  commissioned  value,  which  if  a standard  landing  speed  of  v = 250  ft/sec 
is  assumed,  could  result  in  an  error  on  the  order  of  19  ft/sec . The  obvious  solution  is 
more  frequent  checks  of  the  elevation  angle  and  making  the  most  recent  value  avail- 
able. For  a glide  slope  over  flat  terrain  like  PVO,  the  variation  would  be  due  mostly 
to  seasonal  fluctuations  in  the  water  table  and  several  checks  a year  could  reduce 
the  error  to  less  than  .05°  or  about  4 feet/sec.  Other  sites  are  not  as  fortunate.  A 
4 foot  snow  cover  can  cause  a .4°  change  in  average  elevation  and  an  antenna  used 
at  New  York's  LaGuardia  Airport  designed  specifically  to  minimize  the  effect  of  terrain 
variations  still  exhibits  a .1°  variation  due  to  5 foot  tides?  These  represent  ground 
speed  errors  of  33  and  8 ft/sec  respectively. 


Consider  the  traces  of  glide  path  structure  produced  by  the  Flight  Inspection 
System.  The  test  procedure  involves  visual  position  fixes  at  runway  threshold  and 
at  the  far  end  plus  a vertical  fix  at  threshold  using  the  radio  altimeter.  These 
potentially  contain  sufficient  error  to  cause  the  structure  plot  to  deviate  noticeably 
but  smoothly  from  zero  near  the  time  corresponding  to  runway  threshold  (See  Figure 
5-C).  This  of  course  is  in  addition  to  the  structure  associated  with  the  facility  itself. 
In  the  region  farther  out  substantial  structure  is  apparent  in  the  RGR  and  OKC  plots 
and  a lesser  amount  in  the  PVO  plot.  While  well  within  the  allowed  guidelines 
(±20  micro-amperes)  the  amount  present  at  the  RGR  and  OKC  facilities  is  sufficient  to 
cause  serious  errors  in  the  calculated  ground  speed.  To  see  this,  suppose  that  a pilot 
descends  to  the  - :nway  in  a straight  line  elevated  at  precisely  the  same  angle  as  the 
average  glide  path  and  intersecting  the  runway  at  the  point  abeam  the  glide  slope 
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antenna. 


I 


His  glide  slope  deviation  receiver  should  ideally  indicate  zero,  with  perhaps 
a small  deflection  near  threshold  because  of  the  offset  of  the  antenna.  This 


is  not  the  case  since  the  signal  received  will  be  — 6(/  that  is  the  negative  of  the 
deviation  of  the  actual  glide  path  from  average  alignment.  To  estimate  the  effect 
of  this  on  v^,  the  following  quantity  was  calculated  using  typical  structure  plots 
for  each  of  PVO,  RGR,  and  OKC. 


I 12  tUt)  y(t) 

*2"*1  l J'l  \s.n  ( 0+  6(t)  )J 


The  time  intervals  correspond  to  3500  feet  to  2 nautical  miles  (Region  A)  and  2 
nautical  miles  to  4 nautical  miles  (Region  B).  6#  is  the  derivative  with  respect 

to  time  of  and  its  calculation  involved  a low-order  polynomial  fit  to  6^, 
which  considerably  smoothed  any  of  the  short  term  (1-2  sec)  fluctuations.  The 
quantity  represents  an  "average"  value  of  the  error  induced  in  the  ground  speed 


due  to  glide  slope  "bends". 


ERROR  (ft/sec) 


Region  A 38 

Region  B 43 


As  can  be  seen  from  various  traces,  the  glide  slope  structure  is  usually  oscillatory 
in  nature  with  decreasing  (angular)  amplitude  and  increasing  period  as  distance  from 
the  antenna  increases.  At  distances  on  the  order  of  10,000  feet  from  threshold  the 
contribution  can  be  quite  large  and  relatively  consta 't.  For  Figure  5-B,  the  average 


( 


ground  speed  error  produced  by  the  structure  between  points  I and  II  is  16  ft/sec 
for  a time  interval  of  15  seconds. 

In  a realistic  flight  the  received  glide  slope  signal  ( 6 ) will  be  the  sum  of 
two  parts,  one  representing  the  deviation  of  the  craft  from  the  Ideal  straight  line 
glide  slope  ( 6 ) and  the  other,  as  discussed  above,  representing  the  deviation,  of 
this  idealized  glide  slope  from  the  null  reading  points  of  the  actual  glide  slope  (-£  ' 
(See  Figure  6). 


(10) 


6 = 6-6 

t s 


Glide  slope  deviation  traces  under  turbulent  flight  conditions  will  have  deviations 

from  zero  with  period  comparable  to  or  smaller  than  the  glide  slope  structure's 

period  (See  Figure  4-C).  Thus  an  averaging  or  other  smoothing  algorithm  to  remove 

6 from  the  glide  slope  deviation  signal  will  be  unsatisfactory  under  many  not 

unusual  circumstances.  Unlike  6 and  v , v is  affected  comparatively  little  by 

y x 

turbulence  or  other  excursions  from  the  glide  path  due  mainly  to  the  inertia  of  the 

aircraft  and  its  flight  being  nearly  parallel  (cos  0 = .9986)  to  the  ground.  This 

suggests  attention  be  directed  to  the  possibility  of  smoothing  v , which  is  usually 

relatively  constant.  If  the  ground  speed  is  averaged  or  smoothed  over  a time 

greater  than  the  period  of  the  fluctuations  of  6 , the  associated  error  in  v will 

* x 

diminish.  For  example,  if  the  form  of  the  fluctuation  is  taken  to  be: 


01) 


f (t)  = sin  ^ t 


and  the  smoothing  procedure  is  simply  to  average  over  time  T^,  the  result  will  be: 
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f (f) 

av 


T 

2 rr  T 


2rT 


(1-co* 


In  the  unrealistically  simple  case  when  v^  is  a constant  averaging  will  work 

because  one  can  always  take  T.  large  enough  to  eliminate  the  error  due  to  6 . 

” s 

The  question  becomes:  When  the  true  ground  speed  varies  in  time,  how  much 


of  its  variation  is  lost  for  various  degrees  of  smoothing?  Discussion  of  this  is 
deferred  until  later. 


CHAPTER  IV 


TERRAIN  - INDUCED  ERROR 


Altitude  data  from  the  radar  altimeter  was  usually  available  from  about 


20  seconds  before  crossing  the  runway  threshold  to  the  stop  end  of  the  runway. 


Comparing  it  with  the  inertial  guidance  derived  altitude  and  also,  in  the  case  of 


Atlantic  City,  with  the  theodolite  tracking  data,  it  was  apparent  that  the  difference 


between  the  two  altitudes  was  basically  a slowly  varying  function  of  time,  but 


with  occasional  rapid  changes.  Unfortunately,  United  States  Geological  Survey 


maps  indicate  that  for  RGR,  OKC  and  PVO  the  land  beneath  the  glide  paths  for 


which  radar  altitude  data  is  available  lies  almost  entirely  between  two  contour 


levels  differing  by  ten  feet  and  no  correlation  between  general  topography  and 


radar  altitude  was  possible.  The  sudden  changes  were  found  to  correspond  well 


with  features  on  the  map  such  as  roads  and  creeks.  Part  of  the  traces  in  the 


Atlantic  City  data  was  of  a highly  irregular  nature  and  judging  from  its  relation  to  the 


known  landmarks,  was  measuring  the  altitude  from  high  grass  or  bushes. 


Any  short  period  abberations  in  radar  altimetry  data  can  be  easily  removed 
by  some  form  of  filtering.  The  slower  variations  however,  are  a very  difficult  problem  . 
For  example,  in  much  of  the  RGR  and  OKC  data,  the  inertial  altitude  trace  is  essentially 


a straight  line  indicating  a constant  rate  of  descent: 


h.  = v t 
1 / 


The  radar  data  however,  often  shows  variation  that  conforms  reasonably  well  over 


some  interval  to  the  form: 


v t + A cos 

y 


The  additional  term  is  of  course  the  consequence  of  flying  over  a hill  or  valley. 

In  the  above,  h is  the  distance  the  aircraft  has  descended  below  some  arbitrary  level 
and  t is  the  time  elapsed  since  the  craft  was  at  that  level . V is  the  rate  of  descent 

y 

and  A and  T are  the  amplitude  and  period  of  the  parametrized  undulation.  The 
error  induced  in  the  vertical  velocity  can  be  obtained  by  differentiating  the 
second  term  giving: 


The  values  A = 4 feet  and  T = 10  seconds  are  appropriate  for  the  region  indicated 
on  Figure  4-B.  From  equation  3,  the  consequences  of  this  are  errors  in  v^,  which 
can  be  catagorized  in  3 ways; 


A maximum  instaneous  error  as  great  as  47  ft/sec 
An  RMS  error  of  33  ft/sec 

An  average  error  of  29  ft/sec  over  a time  interval 
as  great  as  5 seconds 


Similar  examples  with  a slightly  smaller  magnitude  but  longer  period  can  be  seen  in 
the  flight  recordings.  The  above  example  was  chosen  because  the  proximity  of  the 
two  altitude  traces  makes  visual  comparison  easier. 


That  this  much  error  occurs  over  land  whose  elevation  falls  between  two  TO 
feet  contours  is  indicative  of  the  difficulties  that  terrain  can  cause.  As  further  proof 
the  survey  map  was  used  to  make  a complete  profile  of  the  land  directly  below  the 
OKC  approach  path  (Figure  7).  Obviously,  there  is  a considerable  variation 


I 


1 


in  the  slope.  To  illustrate,  consider  the  following  four  intervals  for  which  the 
terrain  slope  and  the  error  it  induces  in  the  ground  speed  (assuming  a 250  ft/sec 
approach)  are  given: 


INTERVAL 

LENGTH 

AVERAGE  SLOPE 

SLOPE-INDUCED 

(feet  from  threshold) 

(seconds) 

ERROR  IN  V (ft/se 

0 - 3500  14  -2.9  XI0"*  +1 


3500  - 12500 

36 

-4.2  x io"1: 

+20 

12500  - 16000 

14 

-1 . 1 x io"; 

+54 

16000  - 20000 

16 

+5.0  X 10 

- 2 

This  is  not  a problem  peculiar  to  Oklahoma  City.  Topographic  profiles  of  the  terrain 
under  the  glide  slope  for  the  first  2500  feet  out  from  threshold,  average  slopes,  and 
the  velocity  errors  induced  were  tabulated  for  several  other  major  airports.  To 
illustrate  the  variability  of  the  slope,  the  2500  foot  interval  was  then  divided 
in  half  and  the  calculation  was  repeated  for  both  1250  foot  segments,  each  representing 
about  5 seconds  of  flight  time. 


AIRPORT 


Greater  Pittsburgh 
Dallas  (Love  Field) 
Baltimore  (Friendship) 
St.  Louis  (Lambert) 
Washington  (Dulles) 


AVERAGE  SLOPE/ Vx  ERROR  (ft/sec) 
(2500  feet)  (first  1250  feet)  (last  1250  feet) 


+.004/-  19 
+.002  /-  10 
-.026/ +124 
-.001  / + 5 
+.001  / - 4 


-.  08/ +380 
-.016  / + 76 
-.016  /+  76 
-.010/+  46 
+.007  /-  34 


+.088/ -418 
+.  02  / - 95 
-.036/  +171 
+.007/-  34 
-.006  /+  27 
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Clearly  the  above  represents  a degree  of  error  which  is  intolerable.  As 
with  the  glide  path  irregularities  the  question  of  possible  mathematical  treatment 
of  v^  data  to  supress  the  incorrect  contributions  arises. 

Unfortunately,  the  variation  in  time  of  the  correct  and  incorrect  contributions 
is  often  so  similar  as  to  be  mathematically  indistinguishable.  Examination  of  the 
Arma  inspection  Cata  glide  slope  structure  traces  confirms  the  existence  of  oscillations 
nearly  sinusoidal  in  nature  with  periods  ranging  from  about  30  seconds  of  flight  or 
7500  feet  spatially  down  to  1000  feet  or  4 seconds  (See  Figure  8).  United  States 

• I 

Geological  Survey  maps  exhibit  the  same  variations,  with  an  even  wider  range 

* •. 

in  time  scales  possible.  A simulated  flight  through  a thunderstorm  outflow  produces 
variations  in  both  airspeed  and  ground  speed,  which  are  strikingly  sinusoidal  with 
large  amplitude  (over  30  ft/sec)  and  a period  of  about  20  to  40  seconds  (See  Figure  9).5 
Obviously,  to  devise  an  algorithm  to  remove  the  unwanted  contributions  would  also 
destroy  the  usefulness  of  v^  in  the  above  situation. 

Examination  of  the  terrain  profiles  makes  it  obvious  that  no  simple  average 
slope  of  terrain  correction  will  be  effective,  but  rather  only  the  storage  in  an  on 
board  computer  of  fairly  close  approximation  of  the  actual  profile  will  be  sufficient. 

Even  supposing  this  were  feasible,  there  remain  substantial  difficulties.  Every  landing 
would  involve  taking  a position  fix  several  miles  out  and  then  a continuous  updating 
of  position,  so  as  to  extract  the  correct  terrain  slope  from  the  stored  profile.  Visual 

4 

fixes  are  often  not  available;  using  the  outer  marker  would  involve  several  hundred 
feet  of  error.  Calculating  the  point  of  interception  of  the  glide  slope  using  the  beam's 
geometry  and  the  plane's  altitude  is  unreliable  because  of  uncertainties  in  the 
glide  slope  and  it's  structure  and  lack  of  knowledge  of  the  elevation  above  the  plane 
of  the  glide  slope  antenna.  To  then  use  the  somewhat  inacurate  fix  and  integrate 


the  determination  of  v^  to  update  position  is  an  algorithm  whose  errors  can  increase 
dramatically  under  very  simple  conditions.  In  fact  if  the  terrain  slope  ever 
equalled  the  glide  path  angle  for  some  interval  the  method  would  fail  completely. 

Barometric  altitude  probably  is  more  useful.  Unfortunately,  its  accuracy  and  the 
degree  of  the  problem  associated  with  the  delayed  response  to  changing  descent 
rates  cannot  be  evaluated  since  no  flight  recordings  with  barometric  traces  were  avail- 
able. It  should  be  pointed  out  that  during  typical  wind  shear  conditions  ground 
based  barometers  have  recorded  rates  of  change  in  barometric  pressure  as  great  as 
two  millibars  in  15  minutes?  If  this  is  attributed  to  the  movement  of  a cold  air  outflow 
from  a thunderstorm  past  the  measuring  instrument  at  about  20  knots,  an  aircraft  flying 
at  140  knots  could  encounter  the  same  change  in  1/8  the  time,  or  a rate  of  .018 
millibar/sec,  which  corresponds  to  a ficticious  descent  rate  of  .50  ft/sec,  yielding 
an  average  weather  induced  ground  speed  error  of  about  10  ft/sec  for  a time  interval 
of  over  100  seconds. 
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CHAPTER  V 


CONCLUSION 

Erroneous  contributions  to  the  ground  speed  by  irregularities  in  the  glide 
slope  and  natural  terrain  variations  have  profiles  In  time  so  similar  to  the  actual 
variations  in  ground  speed  that  mathematical  treatment  cannot  effectively  discriminate 
between  them.  The  size  of  the  errors  depends  of  course  on  the  nature  of  the  glide 
slope  and  terrain,  but  from  what  has  been  seer^a  quite  conservative  estimate  is 
that  a typical  approach  will  likely  have  ground  speed  errors  of  at  least  10  ft/sec 
due  to  glide  slope  and  40  ft/sec  due  to  terrain. 
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FIGURE  2 


An  aircraft  at  distance  x from  the  point  O,  and  lying  in  the  plane  tangent  at  O 
to  the  Earth's  surface  is  located  a distance  d above  the  surface: 


d = (R2  + x?)  /2  - R 


For  x very  much  smaller  than  R = 3440  nautical  miles,  the  first  few  terms  of  the  binomial 

\/2 

expansion  of  (R2  + x2)  give  an  excellent  approximation.  For  x and  d in  nautical 


miles: 


d = R(l+  1/2  (Vr)2  +...)_  r*- 
.2 


2R 


= .000145  x2 


For  d in  feet  this  is:  d = .833  xr 

For  x small  compared  to  R the  two  radii  shown  are  very  nearly  parallel,  and  the  correc- 
tion d may  be  added  algebraic ly  to  the  height  above  the  tangent  plane  h to  get  the 
altitude. 
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Figure  4-A  - Anna  Flight-  Recorder  Output  for  ILS-3  Run  for  OKC  on  7-28-77 
The  horizontal  seal*  Is  I second  per  block.  The  vertical  scale  is 
5 ft/bloek  for  altitudes  and  .035  degrees/1) lock  for  the  glide  slope 
deviation  signal.  The  above  scale  is  standard  for  all  further  diagrams. 

The  ground  speed  as  derived  from  the  DME  is  fairly  constant  at  245  ft/sec. 
The  obvious  irregularity  in  spacing  at  point  I is  clearly  physically  impossible. 
Further,  the  vertical  descent  rate  as  derived  from  the  inertial  altitude  is  essentially 
constant,  strongly  suggesting  a near  perfect  descent  with  « = 0.  However,  the 
glide  slope  deviation  signal,  if  true,  implies  a sinusoidal  oscillation  of  the  craft 
with  accellerations  greater  than  2g  occurring  in  Region  II.  Therefore,  the  glide 
slope  deviation  signal  is  in  error,  or  at  least  must  be  measuring  something  other 
than  the  position  of  the  aircraft. 

Note  the  discrepancy  between  radar  and  inertial  altitude.  The  bump 
in  the  radar  altitude  at  about  5 seconds  before  threshold  corresponds  to  a road 
on  topographic  maps  of  the  area,  but  the  tendency  of  the  two  traces  toward  con- 
vergence at  14  seconds  contradicts  the  maps  and  might  be  attributed  to  a change 
in  vegetation  height. 

Figure  4-B  - Arma  Flight  Recorder  Output  for  ILS-3  run  for  RGR  on  7-29-77 
The  variation  in  the  radar  altitude  trace  at  6 seconds  from  threshold 
corresponds  to  a creek  on  the  topographic  maps.  The  valley  is 
referred  to  on  page  13  of  the  text. 
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Figure  4-C  - Arma  Flight  Recorder  Output  for  ILS-3  run  for  PVO  on  11/1 1/77 
This  run  was  made  under  turbulent  conditions  as  indicated  by  the 
variations  in  the  glide  slope  deviation  trace.  Modeling  the  trace  in 
Region  I as  a sine  wave  of  period  12  seconds,  we  can  estimate  that  a 
vertical  accelleration  as  great  as  .2g  was  encountered. 

Figure  4-D  - Arma  Flight  Recorder  Outpu,t  for  ILS-3  run  for  PVO  on  11/11/77 
The  flight  conditions  were  now  calm  as  evidenced  by  the  glide  slope 

and  constant  descent  rate.  Note  however,  the  variations  on  the 

jj 

radar  altimeter.  The  notches  in  the  trace  correspond  well  with 
roads,  but  the  other  features  can  only  be  vegetation  which  does  not 
show  up  on  topographic  maps.  Note  the  uneven  spacing  of  the  DME 
output  in  Region  I . 


| 


27 


FIGURE  5-A 


FIGURE  5-C 
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Figure  5-A  - Arma  Flight  Recorder  Glide  Slope  Structure  ( 6^)  Output  for  OKC 

i 

For  this  and  all  further  plots  of  6$  the  vertical  teale  Is  .05  degree/line 
and  the  horizontal  scale  Is  1 sec/line  or  approximately  245  ft/sec. 

" ; ■ ‘ ’ • ’ . . ■ • j 

Figure  5-B  - Arma  Flight  Recorder  Glide  Slope  Structure  ( 6^)  Output  for  RGR 
Points  I and  It  refer  to  page  10  of  the  text. 

Figure  5-C  - Arma  Flight  Recorder  Glide  Slope  Structure  ( 6 ) Output  for  PVO 

* 

This  is  an  exceptionally  structureless  glide  slope.  The  hypeibolic 

variation  near  threshold  results  from  an  error  in  the  position  fix, 

■ ' 
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Figure  9 - 


Ground  speed  of  a DC -8  with  a tailwind  of  four  (4)  knots  and  controls 
fixed  abruptly  encountering  a headwind  of  about  forty  (40)  knots. 
Vertical  currents  were  also  present.  Initial  qpeed  was  Ve  = 230  ft/sec . 
Note  the  similarity  in  periods  between  these  oscillations  and  the  glide 
’slope  undulations. 
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